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(D ! Abstract 



a 



Radiofrequency vortex spin-transfer oscillators based on magnetic tunnel junctions with very 
low resistance area product were investigated. A high power of excitations has been obtained 



n^ • characterized by a power spectral density containing a very sharp peak at the fundamental frequency 

^, 

O . and a series of harmonics. The observed behaviour is ascribed to the combined effect of spin transfer 

torque and Oersted- Ampere field generated by the large applied dc-current. We furthermore show 

(N 

^ ■ that the synchronization of a vortex oscillation by applying a ac bias current is mostly efficient when 

^O . the external frequency is twice the oscillator fundamental frequency. This result is interpreted in 

OO 

^^ \ terms of a parametric oscillator. 

o. 
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I. INTRODUCTION 

The interaction between a spin-polarised current and a local magnetization, via 
spin transfer torque, can induce steady-state magnetization precessions. This effect is 
intensively studied both from a fundamental point of view as well as for its potential 
application as radiofrequency oscillators. Two classes of these so-called spin torque 
oscillators (STO) are being studied. In a first class, the magnetization of the excited 
layer is close to single domain and precesses around its equilibrium position. In the 
second class, magnetization forms a magnetic vortex-"- that is driven into a periodic 
circular motion by spin transfer torque. In both cases, improving the output power and 
the quality factor of these STO is a challenging requirement for practical applications of 
these devices. In this context, we studied vortex oscillators based on magnetic tunnel 
junctions that have a high output power up to 20 nW. This output power might be 
further enhanced by synchronizing several oscillators. This concept has been tested by 
using either magnetic coupling^ or direct serial or parallel electrical connexion^. It was 
also shown that a STO could be locked on an external source at a frequency close to 
its natural frequency /o, either by feeding the oscillator with a ac bias current^""- or by 
applying a RF magnetic field^. Theoretically, locking should also be possible at any 
frequency being a rational multiple of the natural frequency—. This was experimentally 
demonstrated on macrospin STO^. In this paper we show however that synchronization of 
our vortex oscillator is not obtained for many frequencies. In particular, the most efficient 
synchronization is observed when applying a ac bias current at a frequency equal to twice 
the fundamental frequency. This result is explained by considering the complete system 
formed by the vortex oscillator and the external microwave source as a parametric oscillator. 



II. SAMPLE CHARACTERIZATION 

Our samples are ultra-low RA magnetic tunnel junctions (MTJ) in which vortex 
oscillations can be observed when the junction is subjected to a large dc bias current 
and a low in-plane field. The large dc current produces both a large Oersted-Ampere 
field which, together with the in-plane demagnetizing field, is responsible for the vortex 



nucleation and a spin transfer torque that starts the vortex oscillation. A similar be- 
haviour was observed in low RA MTJ subjected to out-of-plane polarized spin current^. 
In contrast to these experiments, the vortex oscillations are observed in our case 
with in-plane polarized spin current. Our MTJ are composed of the following layers: 
IrMn-j/CoFe2/ RuQ,-j/CoFe2.^/ AlOx/CoFe^/NiFe^ where the subscript corresponds to 
the layer thickness in nm. After deposition, the stack was etched into a pillar with a 
300 nm diameter circular cross section. The MR ratio is about 12% and the resistance 
area product is equal to 0.3 Vt.^m?. In such low RA junctions, there may be some spatial 
inhomogeneities in the current flow due to hot spots in the barrier as suggested in the work 
of Houssameddine et al.— . Four point static resistance measurements were performed on 
the samples as a function of the in-plane magnetic field applied along the easy-axis (R-H 
curves). At low current, the usual sharp transition is observed between the low resistance 
state (parallel state) and the high resistance state (antiparallel state). At higher current 
value (|/| > 5 mA), the R(H) hysteresis loops become more slanted and steps and hysteresis 
appear (Fig. [1]). These features become broader when the bias current is increased. This 
observation is interpreted as the formation of a vortex due to the Oersted-Ampere field 
generated by the current. This assumption is supported by the fact that our R-H curves are 
similar to those observed in thick circular samples where the stable micromagnetic state is 
a vortex—. The vortex state is characterized by an average magnetization close to zero and 
thus by a device resistance equal to the mean value between the resistance values in parallel 
and antiparallel magnetic configuration. By increasing the applied field, we observed first 
a gentle slope of the curve that is ascribed to the vortex core moving toward the edge of 
the sample and then a sharp transition corresponding to the vortex expulsion at large field. 
The larger the applied current, the larger the Oersted- Ampere field and the more stable 
the vortex state. As a result, the low field slope of the R(H) curves decreases and the 
field value corresponding to vortex expulsion increases with bias current. Some hysteresis 
is observed on the curves at high current indicating that besides this overall picture, more 
complex detailed micromagnetic configurations may appear in these pillars. Nevertheless, 
the R-H curves although being dependent on the magnetic preparation of the sample 
before measurement are reproducible provided they are measured in the same conditions. 
The variation of Rmax as a function of current is due to the large bias dependence of the 
resistance in antiparallel configuration. 
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Figure 1: (Color online) Magnetoresistance as a function of the applied field for different currents. 
The curves are vertically shifted in order to overlap in the parallel state. The field is applied in- 
plane along the easy axis. The sample is prepared in the parallel state (—600 Oe) before the field 
cycle. The cycle starts from the parallel state to the antiparallel state and back to the parallel 
state. The electrons fiow from the reference to the free layer. Similar behaviour is obtained for 
positive current. 



III. SPECTRAL MEASUREMENTS 



Spectral measurements were also performed on our samples. The MTJ is polarized 
by a dc current through a bias tee and the output power spectral density is measured 
with a spectrum analyzer. For bias currents larger than 17 7nA (i.e. 2.4 IQ^ A/cm?) in 
absolute value, a radiofrequency spectrum composed of one large and narrow peak at 



the fundamental frequency /o and of about 10 harmonics of much smaller amplitude was 
observed (Fig. [2]). The amplitude of the largest peak is up to 1 fiV'^/Hz in the best 
samples and its linewidth is about 1 MHz which corresponds to an integrated power of 
20nW. The fundamental frequency varies from sample to sample between 350 MHz and 
500 MHz and the relative amplitudes of the various peaks differ from sample to sample. 
The harmonic 2/o is always much larger than the other harmonics and its voltage amplitude 
is about 10 times smaller than the fundamental peak. Its linewidth is twice larger than the 
linewidth of the fundamental peak. In fact, we more generally observed that the linewidth 
of the A^-th harmonic NJq is 2^ times larger than the fundamental peak linewidth. When 
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Figure 2: (Color online) Power spectral density emitted by the junction in which the pre-existing 
static vortex is driven into motion. Measurements are performed with a bias current of —20 mA 
and under 110 Oe in-plane applied field. Each peak was measured separately for a better resolution: 
harmonics 1 to 5. Inset: harmonics 6 and 8 to 12 (the 7th harmonic could not be detected). 



varying the static field or the dc-current, the peak frequency hardly changes showing only 
a slight increase of frequency with increasing the bias current ( lOMHz/mA). This poor 
tunability, compared to niacrospin STO, was already observed in nanocontact-based vortex 
oscillators^ii^ii^. In addition a strong hysteretic behaviour is also observed while sweeping 
the dc-current. When decreasing the bias current (in absolute value), the dynamical vortex 
state survives even at current smaller that the critical current Ion corresponding to the 
onset of the dynamical state. The observation of two threshold currents /„„ and loff that 
correspond respectively to switching on and off the dynamical state was already reported 
in the case of nanocontactsiS. In our samples, Ion and loff are respectively about 17 mA 
and 6 mA (i.e. 2.4 10^ and 8.4 lO^A/cm^). 

The onset of a vortex dynamical state due to spin transfer torque is well understood in the 
case of out-of-plane spin polarized current^. This case has been extensively studied in the 
literature where most experiments are performed with out-of-plane applied G.e\d^^^^^. In 
our experiment however, magnetic field is applied in-plane. In this case, it was pointed out 
that no vortex motion is expected for in-plane spin-polarized current, except in the case of a 
non uniform magnetic distribution in the polarizer—. In our case, the Oersted-Ampere field 
is quite large (250 Oe for I^ic = 20 niA) compared to the RKKY coupling field estimated 
to 500 Oe assuming a reasonable value of the coupling energy {Jrkky = 1 erg/cm'^). 
Therefore the magnetization distribution of the reference layer is probably non uniform, 
being either a C-state or a strongly off-centered vortex. This magnetic distorsion in the 
pinned layer may explain the hysteresis observed at 20mA in Fig. [H This distribution 
is responsible for an inhomogeneously polarized spin current that exerts a torque on the 
vortex core in the free layer leading to a gyrotropic motion of the vortex. When the current 
is reversed, the C-state will be reversed despite the presence of the exchange bias field and 
the current sign will change, so the spin torque term remains rougthly the same. In fact, 
this is experimentally observed: the vortex dynamical state exists for both positive and 
negative current. Finally, let us gain some insight into the vortex trajectory. Considering 
the numerous harmonics observed in the measured spectrum, the vortex orbit may be a 
more complex trajectory than a mere circle, such as an irregular ellipse along which the 
vortex core travels at non-constant velocity. From the total power of the first peak in the 
power spectrum, the resistance change during a period may be estimated and compared 
to the magnetoresistance measured for the corresponding dc-bias current; it is about 10% 

6 



of the maximum possible resistance change. So it is reasonable to assume that the area 
enclosed by the vortex trajectory in the free layer represents 10% of the total surface of the 
sample, which means an average orbit radius of 1/3 of the sample radius. 



IV. SYNCHRONIZATION 

Synchronization of the oscillator with an external signal was then tested by applying 
an additional ac bias current to the magnetic tunnel junction. Since the sample is a one- 
port device, a power-splitter was used in order to connect the sample both to the spectrum 
analyzer and to the microwave source. The sample is therefore simultaneously subjected 
to a dc-current Idc and to a small ac-current lac- When sweeping the frequency f^xt of 
the injected power, we observe that the oscillator fundamental frequency is tuned by the 
external signal. In Fig. [3] we show this phenomenon by sweeping the frequency fext either 
around /q, the oscillator fundamental frequency, or around 2/o. In the latter case, we observe 
that synchronization remains possible even for low injected power. Thus synchronization at 
twice the fundamental frequency appears to be most efficient. In that case, synchronization 
takes place on a range of a few tens of MHz around the oscillator natural frequency which 
corresponds to a relative detuning range A///o of about 4%. This modest detuning range 
may be directly related to the poor tunability of the vortex oscillator-. 

Experimentally, synchronization at the oscillator fundamental frequency is observed on 
the spectrum analyzer as the superposition of the relatively broad peak of the oscillator 
and the very sharp peak of the microwave source. When synchronization is performed 
at twice the fundamental frequency, the oscillator peak is observed alone. Therefore it is 
possible to perform another experiment with keeping the external signal frequency fixed 
within the synchronization range and varying the input power. Fig. H] shows the oscillator 
fundamental peak amplitude as a function of the input power. Above a critical input power 
value (about —25 dBm corresponding approximatively to lac ~ 0.5 mA « Idc)-, the noise is 
drastically reduced around the fundamental peak and the amplitude of this peak increases by 
several order of magnitude. Simultaneously the peak linewidth becomes extremely narrow 
(< IkHz), whereas the integrated power is conserved. 
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Figure 3: (Color online) Normalized frequency of the fundamental peak versus frequency of the 
injected microwave; the power of the microwave is —20 dBm in the upper panel [lac ^mA) and 
-28 dBm in the bottom panel {lac 0.4mA); I^c = -20 mA and H = 250 Oe. 

V. MODEL 



All our experimental observations may be explained by considering the equation of motion 
of the vortex subjected to the spin torque due to both the dc and ac bias currents. In order to 
take into account the time-dependent velocity of the vortex and consequent vortex distorsion, 
the vortex motion is described by the generalized Thiele equation^ including the inertial 



term^: 



MX = GxX - dW/dX - r]X + Fst 
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Figure 4: (Color online) Amplitude of the fundamental peak versus power of the injected microwave 
at 2/0. fext = 700 MHz; I^c = —17 mA; in-plane field H = 15 Oe. The blue line is only a guide 
for the eyes. 

where X is the position of the vortex core, M is the vortex mass, G = —GqZ is the girovector, 
7] is the damping constant and —dW/dX is the force due to the magnetostatic potential 
Since the vortex is confined within the magnetic dot due to the strong Oersted- Ampere field, 
we choose a ad-hoc form of the potentiaL We assume that the potential has a parabolic 
shape in the center of the sample^ and confines the vortex within a limit radius ^: W{r) = 
kC,"^ [2(1 — r"^/^"^)] , so that the corresponding force may be written as —K,f{r)X where 
/(r) = (1 — r^/^^)~^. Finally, the effect of the spin transfer torque is described as a force Fst 
acting on the vortex and proportional to the total applied current J(t) = Idc + lac cos(2a;t) 
(with Uext = 2a;), so that it may be written as Fst = XI{t)z x X, X being a proportionality 
factor—. So the equation can be re-written in the following form: 



MX 

Kf{r) XI dc 
-Mdc nf{r) 
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X = 



It is worth emphasizing that Eq. ([2]) is a non-hnear differential equation, by contrast 
to the worlc of Choi et al.— , even when lac is set to zero, because of the definition of the 
potentiaL Such type of equation is expected for a self-sustained oscillator. This equation 
with lac = may be solved assuming a circular trajectory of the vortex core with a radius R. 
The critical current Ic is defined as the minimum current that compensates the dissipation 
7]. In that case, the oscillator pulsation uq is found to be a real number and the orbit radius 
is close to zero. When the bias current increases, the model predicts that the pulsation cuq 
increases as well as the orbit radius. For sub-critical bias current, there is no steady-state 
motion of the vortex, since the pulsation Uq is a complex number with an imaginary part 
describing the damping of the oscillation. All these predictions are qualitatively consistent 
with experimental observations. Quantitative comparison between model and experiment is 
more difficult due to the fact that some parameter values are not precisely known. Never- 
theless, from the measured resonance frequency it was possible to estimate the vortex mass 
which is found of the order of 10~^^g, consistently with previous theoretical calculation^. 

Now consider the full equation (Eq. [2]) when the ac-current is injected into the oscillator. 
We recognize the equation of a 2D parametric oscillator, similar to a damped Mathieu 
equation. When Uext = Swo (i.e. u = Uo), this equation predicts that a dynamical instability, 
defined as an exponential divergence of the orbit radius, can be reached for a small non-zero 
value of the excitation in the presence of damping. In our model, the confinement potential 
prevents the divergence and a stable periodic motion can be obtained, as experimentally 
observed. Other similar cases with an increasingly larger excitation threshold also exist 
when u = Uq/u with n being an integer i.e. when Uext = 2uj = 2uQ/n. This has also 
been observed experimentally : for Uext = (jJq {n = 2), a larger power is necessary to 
lock the peak and for Uext = i^o/2 {n = 4), no efficient locking could be observed. These 
observations are clearly seen in Fig. [5] which represent the Arnold tongues^. This plot is 
obtained in the following way: for various power of the injected microwave, the locking 



10 



range is measured and the minimum and maximum frequencies of the locking range are 
reported in the graph. The V-shape regions thus obtained (shaded area in Fig. |5]) define 
areas where synchronization is possible. This plot clearly shows that synchronization at 
twice the fundamental frequency is obtained with a much smaller excitation compared to 
synchronization at the fundamental frequency itself. Our model also shows that shifting Uext 
from the 2u!o/n values must be compensated by a larger power: in fact, we observed that 
excitation at Suq requires a larger power and excitations at 3/2ci;o or 4a;o were ineffective 
within the experimental power range. 
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Figure 5: (Color online) Arnold tongues : limits of the locking range of the oscillator as a function 
of the injected power. This graph is constructed from a series of measurements similar to those 
reported in Fig. [3] with varying the power of the excitation. The measurements of Fig. [3] correspond 
to an injected power of 10/iVF and 1.6fxW . Solid lines are guide to the eye. 
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Figure 6: (Color online) Output power as a function of the bias dc-current for different values of 
the ac-current l^c- Currents are normalized to the critical current defined as the minimum current 
value necessary to start the vortex oscillations. For sub-critical currents (shaded area), the model 
predicts that an microwave excitation can switch on the oscillator. 



Finally it is possible to find the stability condition of Eq. ([2]): 

.2 y?i. 



{Mu? + Gow - Kf{R)\ + [r/w - AJ, 



dc\ 



2j2 



(3) 



This formula is obtained in the limit range where the dc-current is much larger than the 
ac-current i.e. when I^c » lac, assuming a circular orbit of the vortex core. From this 
equation, it is possible to estimate the output power of the oscillator (supposed proportional 
to the orbit surface R"^) as a function of the applied dc-current I^c for different values of 
the ac-current lac (Fig. [6]). It is observed that the oscillator output power P rises with 
dc-current due to an increase of spin transfer torque. It is also enhanced by the microwave 
excitation as expected for a parametric oscillator. In particular, for subcritical dc-current. 
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the model predicts that sufficiently large ac-current could start vortex oscillations=ii. 

VI. CONCLUSION 

In this paper we demonstrated that low RA MTJ subjected to low in-plane field and 
high current are powerful vortex oscillators. The vortex oscillation is ascribed to spin 
transfer from a non uniform polarizer onto the vortex core in the free layer. This vortex 
oscillator is particularly efficiently synchronized by an ac current which frequency is twice 
the oscillator natural frequency. All our observations on synchronization can be described 
in the framework of a confined 2D parametric oscillator. Parametric amplification be- 
low the critical current is also expected from the model and has to be checked experimentally. 
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